pa fs 
| | 


人 在 甘 明 二 
第 39 卷 第 9 期 计算 机 应 用 研究 Vl 


月 十 | 
No.9 


录用 定稿 Application Research of Computers Accepted Paper 


基于 FDS 和 元 胞 自动 机 动态 耦合 的 火灾 踊 散 模型 
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(武汉 理工 大 学 a. 中 国 应 急 管理 研究 中 心 ;b. 安全 科学 与 应 和 急 管 理学 院 ， 臣 汉 430070) 


摘 要 : 为 研究 火灾 场景 下 温度 、 烟 气 和 CO 浓度 等 灾害 因子 对 疏散 的 影响 ， 建 立 基 于 FDS 和 元 胞 自动 机 动态 耦合 
的 火灾 足 散 模型 。 模 型 将 FDS 的 网 格 和 元 胞 自动 机 的 元 胞 一 一 对 应 ， 把 由 FDS 运行 得 到 的 灾害 数据 通过 Python 等 
技术 手段 实时 加 载 到 元 胞 中 ， 使 灾害 数据 持续 影响 行人 转移 概率 ， 从 而 实现 灾害 和 足 散 的 动态 耦合 ; 以 单 层 教学 楼 
作为 仿真 场景 进行 模拟 分 析 , 对 火 源 位 置 和 热 释放 速率 等 因素 进行 讨论 , 得 出 这 些 因素 对 行人 疏散 进程 的 影响 规律 ; 
将 模型 与 传统 软件 和 同类 方案 进行 对 比 ， 以 对 比 其 异同 。 研 究 表明 : 火灾 导致 的 高 温和 烟 气 会 影响 行人 对 玖 散 路 径 
和 安全 出 口 的 选择 ; 热 释 放 速 率 越 大 ， 行 人 越 早 处 于 危险 状态 ， 同 时 处 于 危险 状态 的 行人 也 越 多 ; 该 模型 相 比 传统 
足 散 软 件 不 仅 能 考虑 火灾 产生 的 致 灾 因 子 对 行人 疏散 的 动态 影响 ， 还 能 确定 行人 最 早 处 于 危险 状态 的 位 置 和 时 间 ， 
并 用 可 视 化 的 方式 表现 出 来 。 

关键 词 : 火灾 疏散 ; FDS; 元 胞 自动 机 ; 动态 耦合 ; 仿真 模拟 
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Fire evacuation model based on dynamic coupling of FDS and cellular automata 


Li Chao*®?, Li Yufeia。% Huo Feizhoua ?i, Zhang Qinqina。 
(a. China Research Center for Emergency Management, b. School of Safety Science & Emergency Management, Wuhan 
University of Technology, Wuhan Hubei 430070, China) 


Abstract: To study the influence of disaster factors such as temperature, smoke, and CO concentration on evacuation at the 


; fire scene, established a fire evacuation model based on the dynamic coupling of FDS and cellular automata. Mapped the grid 
LO of FDS and the cell of cellular automata one by one, and loaded the disaster data obtained from the operation of FDS into the 
cell in real-time through Python and other technical means, so that the disaster data can continuously affect the pedestrian 
transfer probability, to realize the dynamic coupling between disaster and evacuation; Taking the single-story teaching building 
as the simulation scene for simulation analysis, discussed the factors such as fire source location and heat release rate, and 
obtained the influence law of these factors on pedestrian evacuation process; Compared the model with traditional software 


re and similar schemes to compare their similarities and differences. The research shows that the high temperature and smoke 
caused by fire will affect the choice of pedestrian evacuation path and emergency exit; The greater the heat release rate, the 


earlier pedestrians are in danger, and the more pedestrians are in danger at the same time; Compared with the traditional 


evacuation software, the model can not only consider the dynamic impact of the disaster causing factors caused by fire on 
pedestrian evacuation but also determine the location and time of pedestrians in the first dangerous state, which can be 
displayed visually. 

Key words: fire evacuation; FDS; cellular automata; dynamic coupling; analogue simulation 


0 ”引言 的 应 用 价值 。 但 仿真 软件 也 基于 固定 的 行人 玻 散 模型 ， 用 户 
不 能 改变 疏 散 的 实际 移动 规则 ， 这 是 疏散 软件 最 大 的 局 限 。 
近年 来 ， 建 筑 火 灾 频 发 ， 造 成 了 大 量 的 人 员 伤 亡 和 财产 。 例如, 一 些 学 者 通过 使 用 FDS 计算 出 建筑 发 生火 灾后 的 可 用 
损失 。 据 近 10 年 的 数据 统计 ， 我 国 一 共 发 生 过 3.1 万 起 高 层 ”下 散 时 间 (ASET)， 然 后 再 通过 Pathfinder 等 软件 计算 出 该 建 


建筑 火灾 ， 死 亡 人 数 474 人 ， 直 接 经 济 损失 约 15.6 亿 元 由。 ，” 筑 行 人 的 必须 疏 散 时 间 (RSET), 最 后 通过 比较 RSET 与 ASET 
发 生火 灾后 ， 人 员 在 恐慌 、 人 从众、 冲动 等 心理 特征 下 ， 难 以 的 大 小 来 确定 建筑 是 否 满足 防火 和 琉 散 的 要 求 ， 这 样 就 没 能 


冷静 应 对 形势 并 作出 理性 的 逃离 决策 。 因 此 ， 为 最 大 程度 地 ”很 好 地 考虑 到 火灾 灾害 因子 对 玻 散 过 程 的 影响 扎 习 。 
减少 火灾 造成 的 人 员 伤 亡 ， 开 展 火灾 条 件 下 的 人 员 政 散 动力 鉴于 上 述 原 因 ， 越 来 越 多 的 研究 人 员 将 目光 投放 到 建立 
学 研究 是 非常 必要 的 。 疏散 模型 上 。 总 体 上 讲 ， 玻 散 模型 可 大 致 分 为 宏观 模型 和 微 


考虑 到 火灾 条 件 的 危险 性 ， 在 现实 中 不 可 能 进行 火灾 疏 ”” 观 模型 两 类 中。 宏观 模型 将 人 群 移动 视 为 流体 运动 ， 能 够 高 
散 实 验 ， 因 此 随 着 计算 机 技术 的 发 展 ， 计 算 机 模拟 成 为 一 个 ” 效 计算 出 大 规模 人 群 的 疏 散 时 间 ， 但 宏观 模型 过 于 理想 化 ， 
非常 重要 而 可 行 的 工具 。 国 际 上 开发 的 较为 成 熟 的 行 仆 散 ” ”无 法 反映 个 体 之 间 的 相互 作用 和 异 质 性 。 而 作为 微观 模型 典 
模拟 软件 包括 Pathfinderp 、Building EXODUSBl 和 型 代表 的 元 胞 自动 机 模型 不 仅 能 够 反映 出 行人 个 体 间 的 差异 
FDS+EVACI 等 。 这 些 玻 散 仿真 软件 在 工程 实践 中 具有 较 高 还 能 体现 出 行人 在 朴 散 过 程 中 典型 的 心理 特征 和 行为 反应 ， 
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难以 
将 火灾 数 提 


kK; 


在 一 定 的 时 间 间 


细 三 者 耦合 作用 


上 述 内 容 表 昌 
琉 散 过 程 的 影响 ， 
型 对 火焰 蔓延 和 ; 


实地 再 


中 导入 到 元 胞 


朴 散 效率 。 陈 长 坤 等 器] 为 研究 火 源 对 人 员 玻 散 的 影响 ， 结 合 
元 胞 自动 机 提出 一 种 考虑 多 出 口 吸 引 、 人 员 从 众 行为 与 火 源 
的 场 域 疏 散 模 型 。 
月， 很 多 学 者 已 经 考虑 到 了 火灾 因素 对 行人 
并 建立 了 相应 的 数学 模型 。 然 而 ， 这 些 模 


四 气 扩散 的 描述 过 于 简单 化 和 理想 化 Da， 


14 通 过 FDS 


现实 际 火灾 下 的 玉 G 散 过 程 。Cao 等 
自动 机 中 ， 建 立 了 更 符合 火灾 实际 的 


玻 散 模型 。 但 模型 仍然 存在 场景 设置 比较 简单 (没有 考虑 建筑 


空间 内 障碍 物 结构 对 
程度 较 低 的 问题 。 
机 动态 耦合 上 


1 ”模型 描述 
1.1 


气 扩散 和 行人 玻 散 的 影响 ) 和 


可 视 化 


因此 ,本 文 提出 一 种 基于 FDS 和 元 胞 自动 


的 火灾 玻 散 模型 。 


FDS 和 元 胞 自动 机 简介 


FDS 是 具有 


动 问题 ,可 以 很 好 地 计 


维 可 视 化 动态 功能 的 火灾 模拟 软件 ， 其 基 
于 火 源 驱动 流体 的 CFD 模型 能 够 有 效 描述 低 马赫 数 气体 流 


| 算 火 场 中 的 温度 变化 和 气体 浓度 变化 ， 


既 能 计算 烟 气 流动 和 热 传 递 过 程 , 还 可 用 来 观察 火场 内 的 烟 气 、 


温度 、 能 见 度 、 热 释放 速率 、 燃 烧 产 


掀 浓 度 等 参数 随 火 灾变 化 


而 变化 的 情况 34。FDS 的 基本 思想 是 根据 质量 (组 分 ) 守 恒 、 动 


Op 


+V:pu= 


or 


0 


EPO+V: puut Vp= pg+f +Vr 
t 


@ dp 。 

一 (OH+V:DH = 一 十 dg 一 V.d 十 
可 (名 Pp 4 q+9 
_2RT 
?7 


有 .EO 


: 了 为 气体 密 


度 ， kg /ms; u 为 速度 矢量 ， m/s; 
力 加 速度 ,， m/s*;f 为 外 部 力 矢量 ,N; 


量 守 恒 和 能 量 守 恒定 律 建立 相应 的 基本 方程 ， 其 方程 如 下 : 


(1) 


GO) 


G) 


(4) 
8 为 重 


5 为 牛顿 流体 秋 性 应 


ChinaXiv 合 作 期 刊 


(a) 三 角形 网 格 


Fig. 
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受到 广大 研究 者 的 关注 。Zheng 等 [9 考虑 了 火灾 和 烟 气 对 行 ”后 ， 以 与 疏散 更 新 相同 的 时 间 间 隔 (时 间 步 ) 加 载 到 元 胞 自动 
人 运动 的 影响 以 及 烟 层 场 对 疏散 过 程 的 影响 ， 改 进 了 基于 场 机 模型 中 ， 实 时 变化 的 灾害 数据 通过 动态 改变 元 胞 自动 机 中 
域 模型 的 元 胞 自动 机 模型 。 金 泽 人 等 109 考虑 了 火灾 导致 的 恐 行人 运动 的 综合 场 值 从 而 持续 影响 行人 各 方向 的 运动 转移 概 
慌 心 理 对 行人 移动 方向 的 影响 ， 提 出 基于 火灾 场景 的 元 胞 自 率 ， 这 样 便 实 现 了 每 一 个 时 间 步 灾害 数据 和 踢 散 行为 的 动态 
动机 玻 散 模型 。 萨 木 哈 尔 。 波 拉 提 等 0 使 用 量化 的 趋 近 移 动 耦合 ， 如 图 3 所 示 。 将 灾害 数据 和 下 散 行为 进行 动态 耦合 ， 
强度 来 定义 人 员 移 动 规则 ， 通 过 竞争 占 点 原则 解决 人 员 玻 散 能 够 实时 反映 出 灾害 因子 对 行 仆 散 过 程 的 影响 ， 更 加 符合 
过 程 中 的 冲突 问题 。 江 雨燕 等 (2 通过 分 析 已 有 的 元 胞 自动 机 朴 散 实际 。 

理论 ， 结 合 人 员 玻 散 的 特点 ， 构 建 了 火灾 玻 散 模 型 ， 研 究 表 

明 在 陌生 的 疏散 环境 或 紧急 情况 下 ， 适 当 的 从 众 行为 会 提高 


(b) 矩形 网 格 
常见 的 元 胞 空间 结构 


1 Common cellular spatial structure 


(c) 六 边 形 网 格 


图 1 


有 


J 


泣 


(a) Neumann 型 


(b) Moore 型 (c) 扩展 的 Moore 型 


图 2 常见 的 邻 域 类 型 
Fig.2 Common neighborhood types 
疏散 模型 灾害 模型 
red ss Y 
元 胞 自动 机 FDS 
ES NE 
-对 应 
元 胞 =- | 网 格 
二 十 声 鸭 直到- 生 | 设置 参数 [TF 
个 场 的 权重 系数 ”| 如 | 热 释放 
等 参数 
1 Y 
CAD 图 | Fk 
Watlab 编 程 | | 建立 物理 型 | cm 一。| 建立 物理 模型 | | resin 而 区 
i 
运行 软件 井 对 每 一 
Matlab 编 程 | 确定 行人 运动 规则 个 网 格 的 灾害 数据 
进行 收集 
Y 
将 灭 雷 粗 7 
地租 建立 静态 场 、 动 态 - yy Python 迁 行 到 
Matlab 编 程 | 起 并 用 不 汶 、3 每 个 单位 时 Te 
场 等 基本 玻 散 场 地 所 作 凡 天守 上 据 处 理 
实现 行人 运动 与 火灾 发 展 
在 每 个 时 间 单位 下 的 动态 
耦合 


力 张 量 ，N; 天 为 显 始 ，Jkg;， ?了 为 压力 ， Pa; vv 为 单位 体积 


的 热 释放 速率 , W/m; 
9 为 耗 散 函 数 ; 
对 分 子 质量 。 
元 胞 自动 机 起 源 于 1951 年 冯 。 诺 依 曼 提 出 
理论 ， 在 20 世纪 80 年 代 被 应 用 


2 为 热 通 量 矢 量 , W/m? ; 7 为 温度 ， 


R 为 理想 气体 


常数 ， 丈 为 气体 混合 物 相 


究 中 ， 该 模型 不 同 于 其 


的 物理 方程 或 函数 来 再 
成 。 元 胞 的 状态 是 按照 
更 新 ， 从 而 构成 动态 系统 的 演化 。 常 见 的 元 胞 空间 结 
构 和 邻 域 类 型 分 别 见 图 


行 实时 


1 和 2 所 示 。 


1.2 FDS 和 元 胞 自动 机 的 动态 耦合 模型 
FDS 对 烟 气 流动 和 热 传 递 过 程 的 计算 是 基于 网 格 的 ， 即 


隔 内 ， 可 以 计算 出 每 一 个 


的 生命 游戏 

于 公共 建筑 中 的 行人 玻 散 研 

他 的 动力 学 模型 ， 它 不 是 由 严格 意义 

前 定 ， 而 是 用 一 系列 模型 构造 的 规则 构 
照 构造 的 规则 和 周围 其 他 元 胞 的 状态 进 


网 格 由 火灾 产生 的 


灾害 数据 ， 如 温度 、CO 浓度 、 烟 气 浓度 等 。 同 时 ,元 胞 自动 


等 数据 处 到 


机 朴 散 模型 中 的 行人 运动 也 是 基于 元 胞 (网 格 ) 的 。 
FDS 运行 得 到 的 数据 通过 Python 


EL 


基于 此 , 将 
进行 处 理 


建立 基于 FDS 和 CA 动态 看 
人 的 疏散 模型 
图 3 技术 路 线 图 


Fig.3 Technology roadmap 


值得 注意 的 是 ，FDS 得 到 的 灾害 数据 是 三 维 空间 中 的 ， 
但 元 胞 自动 机 的 行人 运动 空间 为 二 维 平面 ， 因 此 在 设置 的 时 


候 ， 对 数据 进行 了 


这 样 的 处 理 : 


1)FDS 每 个 网 格 的 高 度 为 楼 层 高 度 。 


2) 根 据 建筑 
格 温度 、 烟 气 浓 


生 能 化 防火 设计 的 最 不 利 原则 ， 选 取 每 个 网 


度 和 CO 浓度 的 最 大 值 作为 最 终 导 入 数据 。 


3) 模 型 时 间 步 更 新 频率 设置 为 1 步 /0.25s， 即 每 0.25s 加 


载 一 次 灾害 数据 并 更 新 行人 位 置 。 
基于 FDS 和 元 胞 自动 机 动态 耦合 的 火灾 玻 散 模 型 建立 
网 格 内 , 元 胞 大 小 设置 为 0.4x0.4m( 对 应 地 , FDS 中 的 
四 大 小 也 设置 为 0.4x0.4m), 每 个 元 胞 或 为 空 或 被 墙 


在 二 维 
网 格 平 


瞻 


以 及 行人 占据 。 边 界 为 墙壁 和 安全 出 口 
即 为 疏散 成 功 。 选 取 和 矩形 网 格 作为 朴 散 基 本 单元 ; 选取 


， 人 员 运 动 到 安全 出 
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Moore 型 邻 域 作为 行 仆 散 的 概率 转移 和 矩阵， 以 此 来 确定 行 模型 对 行人 处 于 危险 状态 的 考虑 : 
人 的 运动 方向 。 行 人 的 转移 概率 PB 如 式 (5) 所 示 。 a) 大 多 数 烟 气 中 毒 死亡 的 事故 都 是 由 CO 造成 的 ， 当 行 
B=N"exp(ksSy+tkoDy -KcCy -krTy)1— ny)ey (5) 人 所 在 元 胞 的 CO 浓度 大 于 等 于 500ppm 时 02， 即 可 说 明 行 
N= Dexp(ksSy+koD; -keCy -kT ny)e, (0 人 处 于 危险 状态 。 _ 
> b) 极 端的 温度 也 会 导致 行人 死亡 。 根 据 已 有 研究 4171， 人 
其 中 :， 为 概率 的 正规 化 处 理 ， 确 保 之 妃 =1; ij 分 别 为 目标 在 65C 的 空气 中 无 法 呼吸 ， 因 此 确定 65C 作为 楼 层 火灾 危 
元 胞 的 横 纵 坐标 ; 5; ，D; ,Cy 和 五 分 别 表示 静态 场 、 动 态 场 、 险 临 界 温度 。 在 模型 中 : 当 行 人 目标 元 胞 的 温度 场 五 大 于 等 
烟 气 场 和 温度 场 ，Kks ， 嫩 ，Ke 和 右 分 别 为 5 ，D;，Cy 和 地 于 3.25 时 ， 即 可 说 明 行 人 处 于 危险 状态 。 
的 权重 敏感 系数 ， 反 映 了 不 同 场 强 对 总 场 强 的 贡献 大 小 ， 且 本 文 £ ，k，kec 和 刁 的 取 值 如 表 1 所 示 。 
有 名 ，Kkp，kc ，kr e[0,%]; 加 和 分别 表征 了 其 余 行人 和 障 表 1 本 文 的 参数 取 值 
碍 物 的 占据 信息 ， 其 取 值 方法 分 别 如 式 (7) 和 (8) 所 示 。 Tab.1 Parameter value in this article 
,1 目标 元 胞 被 墙 或 障碍 物 占据 (7) 参数 。” 取 值 说 明 
””【0， 目 标 元 胞 无 墙 或 障碍 物 占据 ks 1 静态 场 为 行人 运动 基本 场 
sh 标 元 胞 被 行人 占据 (8) ko 0.01 行人 越 熟悉 环境 ， 动 态 场 值 应 越 低 
” 【0， 目 标 元 胞 没 被 行人 占据 ke 10000 使 烟 气 场 值 的 数量 级 与 静态 场 一 臻 
静态 场 % 表示 出 口 对 行人 的 吸引 作用 ， 知 疏散 场景 中 不 kr 1 温度 场 值 已 作 处 理 ， 系 数 取 1 即 可 
存在 障碍 物 ， 其 大 小 直接 由 距离 式 (9) 计 算 即 可 ,但 如 果 场 景 1.3 演化 更 新 规则 
中 存在 复杂 障碍 物 ( 如 “U” 型 障碍 物 ) 时 ， 直 接 使 用 距离 公式 模型 采用 并 行 更 新 的 规则 对 场景 中 所 有 行人 的 位 置 进行 
可 能 会 出 现行 人 因 陷 入 局 部 最 优 的 陷阱 而 无 法 行动 的 状况 。 更 新 ， 具 体 更 新 规则 如 下 : 
解决 此 问题 的 常用 方法 是 采用 Dijkstra 算法 计算 各 元 胞 位 置 a) 初始 化 行人 分 布 。 根据 膨 加 场 强 ， 运用 相应 公式 求 出 
与 出 口 的 距离 ， 本 文 设置 的 场景 存在 障碍 物 ， 故 采用 该 算法 行人 转移 概率 ， 并 确定 行人 下 一 步 的 位 置 。 
对 静态 场 进行 求解 。 b) 当 多 个 行人 竞争 同 1 个 元 胞 时 ， 随 机 等 概率 地 选择 1 


5 = max{ min VO Dr 0 7 
yD) (dia) 


上 位 行人 进入 元 胞 ， 竞 争 失败 的 行人 留 在 原 地 。 
(9) c) 确定 这 一 时 间 步 处 于 危险 状态 的 行人 ， 并 将 其 标记 。 
d) 更 新 下 一 时 间 步 所 有 行人 位 置 , 并 将 位 于 出 口 位 置 的 


min VO +0 


其 中 ， 必 , 产 ) 为 不 同 出 口 位 置 坐标 ， 其 中 大 是 出 口 数量 ，e 表 行人 从 场景 中 移 除 。 
示 第 个 出 口 。 e) 重复 步 又 一 e) 直 到 所 有 行人 朴 散 完毕 。 
动态 场 D; 则 通过 参考 其 他 行人 走 过 的 线路 信息 ,描述 行 es 

人 间 的 从 众 行为 ， 表 达 式 如 式 (10) 所 示 。 2 “仿真 模 拟 及 结果 分 析 
Dy; =((1—dif).(1-dec). D1 (10) 2.1 朴 散 系统 的 物理 环境 假设 
《dd-dec)18) .sumDP + 中 一 中 琉 散 仿真 场景 设置 为 48.4X15.2m 的 单 层 教学 楼 ， 对 应 
sumD? = DE + DH + DE + GD 的 FDS 网 格 数 和 元 胞 空间 大 小 均 为 121X86 网 格 (元 胞 )， 如 
DF + Dri + DF) + Dj + Dj 图 4(a) 和 4(b) 所 示 ， 教 学 楼 层 高 3.6m。 设 置 两 个 分 别 位 于 场 


其 中 ， dif 和 dec 分 别 表示 扩散 和 衰减 系数 ， 这 里 均 设置 为 景 左右 两 侧 的 安全 出 口 ， 每 个 安全 出 口 宽度 为 2.4m。 教 室 中 
0.309。 在 初始 时 间 ， 所 有 元 胞 的 动态 场 值 为 0。 每 当 有 人 通 存在 大 量 桌 椅 , 发 生火 灾 时 热 释 放 速 率 (HRR) 符 合 of 增长 规 
过 元 胞 (让 时 ，Dy=Dj+1。d, 和 dd; 是 两 个 修正 系数 ， 当 上 一 律 甸 0， 依据 最 不 利 原则 设置 火灾 增长 类 型 为 快速 火 ， 火 灾 增 
个 时 间 步 元 胞 G7) 为 空 且 当 前 时 间 步 存在 行人 时 ，41=1, 否 长 系数 取 0.04689，HRR 设置 为 5000kw/m2， 初 始 火 源 面积 
则 a=0; 当 上 一 个 时 间 步 元 胞 (ij) 存在 行人 且 当 前 时 间 步 也 为 0.16m2， 火 源 位 置 详 见 图 4(a)。 在 FDS 中 桌 椅 材 质 设 为 
存在 行人 时 ，4;=1， 否 则 d=0; “wood pine ”， 燃 烧 反 应 设置 为 “wood, the yield of soot 
在 真实 的 火灾 场景 烟 气 的 扩散 比 火 焰 的 蔓延 更 快 、 ys=0.015g/g and the yield of CO ys=0.004g/g”， 设 定 墙 体 、 地 
更 严重 。 一 般 情 况 下 ， 火 灾 导 致 的 其 他 威胁 因素 如 结构 失效 板 和 楼 板 为 不 可 燃烧 的 惰性 材料 。 元 胞 自动 机 场景 的 黑色 元 
等 ， 不 会 先 于 烟 气 达到 危险 状态 。 据 统计 ， 在 建筑 火灾 中 约 胞 代表 桌 椅 和 墙壁 ， 绿 色 元 胞 代表 固定 分 布 的 行人 ， 数 量 为 
75% 一 85% 的 死亡 是 由 烟 气 导致 的 801。 烟 气 场 2 表示 烟 气 浓 720。 为 了 减 小 随机 误差 ， 所 有 朴 散 时 间 步 数 均 取 20 次 模拟 


度 对 行人 的 排斥 作用 , 其 值 直接 由 FDS 导出 。 值 得 说 明 的 是 ， 的 平均 值 。 
烟 气 浓度 与 能 见 度 成 反比 例 关 系 ， 如 式 (12)08] 所 示 ， 使 用 烟 在 Pathfinder 中 选用 “Steering” 模 式 ， 设 置 行人 的 肩 宽 
气 场 即 可 表征 烟 气 浓度 对 行人 琉 散 的 影响 ， 故 本 模型 不 再 引 为 0.4m, 行人 行为 设置 为 “Go to any exit”, 最 大 速度 为 1.6m/s。 
入 文献 [1 和 提出 的 能 见 度 场 。 Pathfinder 的 场景 图 如 图 4(c) 所 示 。 
R=C/(M,:K,)=C/K (12) ”2.2 灾害 因子 对 下 散 过 程 的 影响 

其 中 ，R 为 能 见 度 ; C 为 经 验 常 数 ，M, 为 烟 气 粒子 的 质量 浓 120 时 间 步 下 的 疏散 过 程 和 烟 气 扩散 过 程 如 图 5 所 示 ， 
度 ; ,为 烟 气 粒子 的 消光 率 。 不 考虑 火灾 时 ， 教 学 楼 内 的 疏散 人 员 以 最 短路 径 向 教学 楼 两 

火灾 造成 的 高 温 也 会 对 行人 玻 散 造 成 严重 威胁 ， 行 人 在 ，” 侧 的 安全 出 口 进 行 疏 散 , 如 图 5(a) 所 示 ; 当场 景 发 生火 灾 时 ， 
疏散 过 程 中 会 基于 本 能 远离 高 温 区 域 ， 向 温度 较 低 的 区 域 玖 。”” 左 侧 出 口 处 的 温度 和 烟 气 浓度 较 高 , 其 烟 气 扩散 状态 如 图 5(b) 


散 ， 因 此 用 温度 场 来 表示 温度 对 行人 的 排斥 作用 ， 具体 计 ”所 示 ， 此 时 更 多 的 行人 分 布 在 走廊 右 侧 ， 如 图 5(c) 所 示 ， 其 
算 方式 如 式 (13) 所 示 。 中 灰色 元 胞 代表 烟 气 ， 颜 色 越 深 ， 浓 度 越 大 。 

n= (13) 高 温和 烟 气 影响 了 场景 内 行人 的 出 口 选择 。 相 比 图 5(a)， 

1 图 5(c) 中 的 行人 因 左 侧 出 口 温 度 和 烟 气 浓度 较 高 ， 更 多 地 选 


中 ，7 为 温度 数据 ， 由 FDS 导出 ， 并 由 Python 导入 到 元 择 向 右 侧 出 口 进行 疏散 ， 这 导致 右 侧 出 口 严重 拥堵 ， 也 因此 
胞 自动 机 中 ; 友 为 环境 温度 ， 取 20C 。 增加 了 总 体 的 疏散 时 间 。 
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间 步 后 就 无 人 通过 了 
2.4 ” 热 释 放 速 率 对 和 于 人 处 于 危险 状态 人 影响 

为 了 探究 不 同 热 释放 速率 对 行人 处 于 危险 状态 的 影响 ， 
在 场景 1 设置 三 种 不 同 的 热 释 放 速 率 进行 对 比 研究 : (a) 
HRR=500kw/m?, (b) HRR=2000kw/m?, (c)HRR=5000kw/m?。 
不 同 HRR 下 处 于 危险 状态 的 最 大 人 数 和 最 早 时 间 如 图 8 所 示 。 


i 人 


(b) 元 胞 自动 机 场景 网 格 图 
"Hitt | 


本 本 EON WD Een [EE 


疏 散人 数 /人 
全 
3 


ra dl 


EE rr a 

EE 3 
EF nn [ET 
(c)Pathfinder 场景 图 100 上 


图 4 仿真 场景 图 en 
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Fig.4 Simulation scenario 芒 散 时 间 步 / 步 

图 7 不 同 场景 下 玻 散人 数 随 玻 散 时 间 步 的 变化 图 


Fig.7 Variation of evacuation number and evacuation time steps in 


全 全 人 


different scenarios 


(a) 无 火灾 情景 下 的 120 时 间 步 时 的 朴 散 过 程 图 


上 一 一 -- 


时 间 / 秒 


HRR=500 HRR=2500 HRR=5000 


图 8 不 同 HRR 下 处 于 危险 状态 的 最 大 人 数 和 最 早 时 间 
Fig.8 Time and number of people in the first dangerous state and the 
first time under different HRR 


(oO) 火灾 情景 下 的 120 时 间 步 时 的 朴 散 过 程 图 从 图 8 可 以 看 出 ， 热 释放 速率 越 高 ， 行 人 就 越 早 处 于 危 
图 5 120 时 间 步 时 的 朴 散 过 程 及 烟 气 扩散 图 险 状 态 ， 同 时 处 于 危险 状态 的 行人 也 越 多 。 这 是 因为 热 释 放 
Fig.5 Evacuation process and smoke diffusion diagram at 30 time steps 速率 升 高 导致 场 景 内 温度 升 高 速率 变 快 、 烟 气 产量 变 大 ， 更 
2.3 火 源 位 置 对 玻 散 的 影响 多 靠近 左 侧 出 口 的 上 下 两 个 教室 中 的 行人 因 教 室 门 宽度 限制 
为 了 探究 火 源 位 置 对 朴 散 的 影响 ， 现 分 场景 进行 比较 分 ”而 无 法 及 时 玻 散 。 
析 。 每 个 场景 对 应 一 个 火 源 位 置 ， 其 他 参数 均 相 同 ， 如 图 6 除 此 之 外 ， 模 型 还 可 将 某 一 时 间 步 下 处 于 危险 状态 的 行 
所 示 。 人 用 可 视 化 的 方式 直观 地 表现 出 来 ， 9 为 场景 3 在 
i T HRR=5000kw 下 40 时 间 步 时 处 于 危险 状态 的 行人 标记 图 ， 


其 中 红色 元 胞 代表 处 于 危险 状态 的 行人 。 


图 6 不 同 场景 下 的 火 源 位 置 分 布 图 


Fig.6 Fire source location distribution in different scenarios 


图 7 显示 了 3 个 不 同 场景 下 左右 出 口 的 疏散 人 数 随 时 间 图 9 40 时 间 步 处 于 危险 状态 的 行人 标记 图 
的 变化 规律 。 由 图 7 可 知 , 场景 3 的 总 玉 散 时 间 最 少 , 160 时 Fig.9 Pedestrian marker diagram with a dangerous 40 time steps 
间 步 时 行人 就 基本 下 散 完成 ;其 次 是 场景 2， 行 人 完成 疏散 ” 2.5 与 传统 软件 和 同类 方案 的 对 比分 析 
需要 200 时 间 步 ; 场景 1 的 总 下 散 时 间 最 长 ，260 时 间 步 时 Pathfinder 的 模拟 过 程 如 图 10 所 示 ， 与 图 5(a) 相 比 ， 在 
行人 才 全 部 疏散 完成 。 这 主要 是 由 不 同 场景 下 两 侧 出 口 的 利 同样 的 时 间 步 下 , 其 场景 内 存在 更 多 行人 。Pathfinder 与 本 文 
用 率 不 同 导致 的 ， 当 火灾 发 生 在 中 间 位 置 时 ， 行 人 向 两 侧 疏 ”模型 (无 火灾 时 ) 疏 散人 数 随时 间 的 变化 对 比如 图 11 所 示 。 本 


散 , 两 个 出 口 都 能 都 得 到 较 好 的 利用 (图 7 中 场景 三 左 出 口 和  ” 文 的 元 胞 自动 机 模型 使 用 了 Dijkstra 算法 计算 静态 场 ， 

右 出 口 两 条 曲线 相差 不 大 ); 火 源 位 置 越 靠 近 左 侧 的 出 口 , 左 ”Pathfinder 的 路 径 规划 基于 A* 算 法 PU， 而 A* 算 法 是 Dijkstra 
侧 的 出 口 就 越 不 安全 ， 越 多 的 行人 会 选择 向 右 侧 出 口 疏 散 ， 算法 的 扩展 C， 在 本 质 上 二 者 相差 不 大 。 因 此 ， 造 成 差异 的 
从 而 导致 了 右 侧 出 口 的 拥堵 ， 例 如 ， 场 景 一 左 出 口 在 120 时 ， 可 能 原因 是 :Pathfinder 相 比 元 胞 自动 机 更 多 地 考虑 了 行人 之 


202205.00046v1 


chinaXiv 


间 的 碰撞 和 速度 的 变化 ， 如 图 12 所 示 。 

然而 , Pathfinder 最 大 的 局 限 在 于 无 法 考虑 火灾 对 朴 散 过 
程 的 动态 影响 , 在 场景 一 HRR=5000kw/m? 的 情况 下 , 若 用 传 
统 软 件 对 该 建筑 进行 人 员 玻 散 安 全 性 的 评估 : 其 RSET 为 
97.3s，ASET 为 56.7s，RSET>ASET， 行 人 不 能 安全 下 散 ， 
如 表 2 所 示 。 


He yy 


HT 


图 10 ”120 时间 步 时 的 Pathfinder 模拟 过 程 图 
Fig. 10 Pathfinder simulation process at 120 time steps 


一 0 一 元 胞 自动 机 模拟 结果 
一 9 一 Pathfinder 模 拟 结果 
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疏散 时 间 步 数 / 步 


图 11 Pathfinder 与 本 文 模型 疏散 人 数 随 时 间 的 变化 对 比 
Fig. 11 Comparison of changes of evacuation number with time 
between Pathfinder and this model 


但 该 评估 方法 只 是 将 两 种 软件 分 别 运 行 ， 对 所 得 到 的 数 
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据 进 行 比 较 ， 并 没有 很 好 考虑 到 火灾 在 发 展 过 程 中 对 人 员 琉 
散 路 径 和 政 散 行为 的 动态 影响 。 此 外 , 虽然 文献 [12] 也 构建 了 
基于 元 胞 自动 机 的 火灾 朴 散 模型 并 考虑 了 火 源 对 行人 的 排斥 
作用 , 但 该 模型 对 火焰 蔓延 的 设置 过 于 理想 化 (Moore 邻 域 随 
机 扩散 ) 且 没有 考虑 到 烟 气 扩散 和 温度 变化 对 行人 玻 散 的 影 


响 。 


在 本 文 提出 的 基于 FDS 和 元 胞 自动 机 动态 耦合 的 火灾 玻 


= 


散 模型 中 ， 火 灾 造 成 的 灾害 因子 通过 影响 综合 场 值 从 而 影响 
了 行人 的 路 径 和 出 口 选择 行为 ， 该 模型 不 仅 能 够 判断 行人 是 


人 | 
否 自 


安全 下 散 ， 同 时 还 能 够 确定 行人 最 早 处 于 危险 状态 的 位 


置 和 时 间 , 并 用 可 视 化 的 方式 表现 出 来 。 本 文 模型 与 文献 [12] 


模型 的 具体 对 比如 表 3 所 示 。 
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图 12 120 时间 步 时 的 Pathfinder 行人 速度 分 布 图 
ig. 12 Pathfinder pedestrian speed distribution at 120 time steps 
表 2 传统 软件 的 疏散 安全 性 评估 结果 


Tab.2 Evacuation safety evaluation results of traditional software 


灾害 因子 标准 致 灾 时 间 ASET RSET 安全 性 
走廊 能 见 度 <Sm 56.7s 
气 层 高 度 <1.5m - 
出 口 温度 >65C 30.9s 56.7s 97.3s ”危险 
右 出 口 温度 >65°C - 


走廊 CO 浓度 ”>500ppm = 


表 3 本 文 模型 与 文献 [12] 模 型 的 具体 对 比 


Tab.3 The model comparison between this paper and literature [12] 


对 比 项 目 本 文 模型 文献 [12] 模 型 
对 出 口 吸 引 的 考虑 引入 基于 Dijkstra 算法 的 静态 场 引入 基于 Dijkstra 算法 的 出 口 吸引 力 
对 行人 从 众 的 考虑 参考 Kirchner A 等 人 19 的 经 典 动态 场 根据 视野 范围 各 方向 行人 数量 进行 计算 
对 火 源 排斥 的 考虑 综合 反映 为 温度 和 烟 气 浓度 对 行人 的 排斥 引入 基于 欧 氏 距离 的 火灾 排斥 力 
对 温度 排斥 的 考虑 导入 FDS 数据 ， 引 入 温度 场 = 


BR 5 


对 烟 气 排斥 的 考 导入 FDS 数据 ， 引 入 烟 气 场 
对 危险 状态 的 考虑 设置 高 温和 高 CO 浓度 会 使 行人 处 于 危险 状态 的 规则 5 
可 视 化 程度 实现 了 火焰 蔓延 、 烟 气 扩散 和 行人 处 于 危险 状态 的 可 视 化 实现 了 火焰 蔓延 的 可 视 化 


3 ”结束 语 


本 文 为 研究 火灾 下 的 行人 玻 散 提供 了 一 个 新 的 方法 : 建 
立 基 于 FDS 和 元 胞 自动 机 动态 耦合 的 火灾 玻 散 模型 后 , 可 
用 于 不 同 的 建筑 场景 ， 通 过 改变 火 源 位 置 、 燃 烧 材 料 、 热 
放 速 率 和 行人 密度 等 参数 ， 找 到 最 不 利 点 下 的 行人 处 于 和 危 
状态 的 时 间 ， 以 评估 发 生火 灾后 该 建筑 中 的 行人 能 否 安全 
散 ， 为 建筑 性 能 化 防火 设计 提供 参考 。 研 究 表 明 : 

a) 火 灾 导 致 的 高 温和 烟 气 会 影响 行人 对 出 口 的 选择 ， 随 
着 火灾 的 发 展 ， 相 比 更 短 的 距离 ， 行 人 更 愿意 选择 温度 和 烟 
气 浓度 较 低 的 路 径 和 出 口 。 

b) 热 释放 速率 越 大 ， 行 人 就 会 越 早 处 于 危险 状态 ， 同 时 
处 于 危险 状态 的 行人 也 越 多 。 

c) 本 文 提出 的 模型 相 比 Pathfinder 能 考虑 火灾 产生 的 致 
灾 因 子 对 行 仆 散 的 动态 影响 ， 除 了 能 评估 行人 是 否 能 安全 
玉 散 外 ， 还 能 准确 得 出 行人 处 于 危险 状态 的 位 置 和 时 间 。 

今后 在 本 文 提 出 的 耦合 模型 的 基础 上 既 可 以 加 入 前 人 提 
到 的 铀 惫 前 进 等 行为 ， 还 可 以 考虑 因 火灾 导致 的 行人 恐慌 心 
里 。 此 外 ， 可 以 进一步 考虑 多 高 层 建筑 火灾 琉 散 的 建 模 与 分 


Ba 


TH 


析 ， 同 时 也 可 以 将 行人 的 预 动作 时 间 以 及 建筑 的 自动 喷 水 灭 
火 系统 等 消防 性 能 也 考虑 在 模型 中 。 
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